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Weakly nonlinear stability of regimes of free hydromagnetic thermal convection in 
a rotating horizontal layer with free electrically conducting boundaries is considered 
in the Boussinesq approximation. Perturbations are supposed to involve large spatial 
and temporal scales. Applying methods for homogenisation of parabolic equations, 
we derive the system of amplitude equations governing the evolution of perturbations 
under the assumption that the a— effect is insignificant in the leading order. The 
amplitude equations involve the operators of anisotropic combined eddy diffusivity 
correction and advection. The system is qualitatively different from the system of 
mean-field equations for large-scale perturbations of forced convective hydromagnetic 
regimes. It is mixed: equations for the mean magnetic perturbation are evolutionary, 
all the rest involve neither time derivatives, nor the molecular diffusivity operator. 
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1. Introduction 

This paper is a continuation of Zheligovsky (2008), further referred to as Zh08, where 
the relevance of the problem considered here for geophysical applications is discussed 
and the work of other authors in this area is reviewed in detail. 

Our ultimate motivation is the same as in Zh08: we would like to find exam- 
ples of "real" convective dynamos, whose weakly non-linear perturbations have upon 
saturation a complex structure and exhibit a complex temporal behaviour, but re- 
main "stable" being uniformly bounded in time. In such examples various "eddy" 
effects, representing the integral influence of the small-scale dynamics on the large- 
scale structures, on the one hand, should play an important role in the evolution of 
the CHM system; on the other, they must not be introduced ad hoc empirically, but 
rather follow from the basic principles. Saturation is impossible if eddy diffusivity is 
negative. Thus, although we use the same mathematical approach, our aim is in a 
sense opposite to that of Dubrulle and Frisch (1991), Zheligovsky et al. (2001) and 
Zheligovsky (2003), who investigated the effects of negative eddy viscosity and dif- 
fusivity as the mechanisms of hydrodynamic and magnetohydro dynamic instability. 
The perturbations can only have a complex spatial structure and/or demonstrate a 
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complex temporal behaviour, if the amplitude equations that are governing them are 
themselves complex enough; this is one of the reasons, why we are obliged to consider 
a rather general case. 

Weakly nonlinear stability to large-scale perturbations of a free convective hyd- 
romagnetic (CHM) regime V, H, T in a rotating horizontal layer is considered here. 
Our immediate goal is to derive amplitude equations governing the evolution of the 
perturbations under the assumption that the a— effect is insignificant in the leading 
order. The amplitude equations (47), (48) and (38) (the latter equation emerges 
if the perturbed CHM regime is non-steady), that we derive in this paper, involve 
a linear operator of combined eddy dijfusivity correction and a nonlinear operator of 
anisotropic combined eddy advection. Both operators are anisotropic. The system has 
new qualitative features, compared to the system of mean-field equations for large- 
scale perturbations of forced CHM regimes derived in Zh08. It does not involve mean 
flow perturbations, and it is mixed: whereas equations (48) for the mean magnetic 
perturbation are evolutionary, the remaining ones, (47) and (38), do not involve time 
derivatives or molecular diffusivity operators. The derivations follow the way paved 
in ZhOS; consequently, the two papers are similar in structure. 

2. Equations of hydromagnetic thermal convection, 
boundary conditions and the Unearisation operators 

In the Boussinesq approximation the evolution of the perturbation satisfies the vor- 
ticity 

dft 

— = uV^n + V X (V X ^^ - H X (V X H)) + V X (V X res + /JTeg), (la) 
ot 

magnetic induction 

— =77V'H + V X (V X H), (lb) 

and heat transfer 

^ = kv't - (V • v)r 

equations, as well as the solenoidality conditions 

V • V = 0, (Ic) 

V-H = 0. (Id) 

(Following Zh08 we use the vorticity equation in place of the momentum equation to 
eliminate pressure.) The following notation is used: V(x, t) and 

f2 = V X V (le) 

are velocity and vorticity, respectively, of a fiow of an electrically conducting fluid, 
H(x, t) magnetic fleld, T(x, t) temperature, t time, u, rj and k, kinematic, magnetic 
and thermal molecular diffusivities, respectively, r/2 angular velocity of the rotation, 
/9Te3 the buoyancy force, is the unit vector along the axis of the Cartesian 
coordinate system corotating with the layer. 



2 



The CHM regime is free in the sense that no source terms are present in the 
equations (1); consequently, the regime is spatially and temporally invariant, i.e. 
there is no dependence of the coefficients in these equations neither on the position 
in space, nor on time. 

The so called small-angle instabihty problem (see Podvigina 2008a, 2009 and 
references therein) is a well-known example of a problem of stabihty to large-scale 
perturbations, in which linear stability of amagnctic convcctive regimes (rolls) near 
the onset is considered. Note, however, that it does not fall into the class of problems 
considered here (even if the difference between linear and weakly nonlinear kinds of 
stabihty is neglected), since it involves two different large scales; like in Zh08, we 
consider only one large scale. 

The following conditions are assumed on the horizontal boundaries of the layer: 
• no stress: 



dVi 



X3=±L/2 



X3=±L/2 



0, 



X3=±L/2 



xz=±L/2 



X3=±L/2 



= 0, 



perfectly electrically conducting boundaries: 



dHi 



dHo 



X3=±L/2 

isothermal boundaries: 

T 



= 0, 







0; 



X3=±L/2 



X3=-L/2 



X3=±L/2 



X3=L/2 



X3=±L/2 



To. 



= 0; 



(2a) 



(2b) 



(2c) 



(2d) 



(Cartesian vector components are enumerated by the subscript.) 

It is convenient to introduce a new variable = T — 7^ + 5{x^ + L/2), which 
satisfies the equation 

_ = «:V2e-(V-V)e + 5V3 (3) 



and homogeneous boundary conditions: 

e 



0. 



(4) 



X3=±L/2 

We define the spatial mean and the oscillatory part of a scalar or vector field /: 

rL/2 ft/2 



1 1-1/2 fi/2 fl/2 

- / / /(x)da:ida:2dx3, {/}^ /-(/), 

£^ J-L/2 J -1/2 J -1/2 



(5a) 



(/) = lim , 

l^QO J-L/2 J -1/2 J -1/2 

and the spatio-temporal mean and oscillatory parts of / are 

((/)) ^ hm i A/(x, t)) dt, m^f- if); (5b) 

t->oo t Jo 

{f)k and denote the k—th Cartesian components of the respective means. 
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(f)„ = (/)3e3, {f}, = f-(f),; 



)3e3, m.^i- 
(f),= (/)iei + (/)2e2, {f}, = f-(f),; 
((fK ^ ((/))iei + ifhe,, mn ^ f - ((fk- 

Following Zh08, the subscripts v and h are thus used to denote vertical and horizontal 

components of three-dimensional mean vector fields. The superscripts v and h have 
a different meaning: they denote the flow velocity and magnetic field components of 
10-dimensional vector fields {uj, v, h, 9), and 7-dimensional (w, h, 9). 

Linearisation of equations (la), (lb) and (3) in the vicinity of the CHM regime 
V, H, © gives rise to the operators 

Ciuj, V, h, 9) = — — + vV^uj + Vx (Vxu? + vxri 

at ^ 

-H X (V X h) - h X (V X H)) + T— + p\/9 X 63, (6a) 
dh 

C^iv, h) = -— + TyV^h + Vx(vxH + Vxh), (6b) 
at 

39 

£^(v, 9) = -— + kV^9 - (V • V)^ - (v • V)e + 5v3. (6c) 

The operator of linearisation of the system (1) is £ = {C^, C'^, C^). 

We consider weakly nonlinear regimes of perturbations, the amplitude of which is 
of the order of e. The perturbed state O + eu?, V + ev, H + eh, Q + £9 also satisfies 
(1) and (3), and hence the profiles of perturbations a;,v,h, (henceforth called 
perturbations) satisfy the equations 

£^(a;, V, h, ^) + eV X (v X a; - h X (V X h)) = 0, (7a) 

£'*(v, h) + £V X (v X h) = 0, (7b) 

£^(v,^) -£(v V)^ = 0, (7c) 

V • h = 0, (7d) 

V • u; = V • V = 0, (7e) 

V X V = w. (7f) 

Averaging of the horizontal component of the Navier-Stokes equation over the fluid 
volume yields 

^ = (v). X re3 - {Vp)n. 

We assume that there is no mean horizontal fluid flow through the layer, i.e. at any 
time (Vp)^ = and 

(v). = 0. (7g) 
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3. Lcirge-scale wecikly nonlinear perturbations of CHM regimes 

We introduce the slow spatial, X = £{xi,X2), and temporal, T = eH, variables. (We 
assume that the layer width L is of the order of unity, and hence a slow variable in 
the vertical direction is not introduced.) For vector fields depending both on fast and 
slow variables, integration in the definitions (5) of mean fields is performed only over 
the small-scale fast variables x or (x, t). The exponent in the temporal scale ratio is 
suitable for CHM regimes with an insignificant a— effect. A power series solution to 
the problem (7) is sought: 



oo 



u:^J2 '^"(x, t, X, T) v^Yl Vn(x, t, X, T) e^, 

n=0 n=0 



oo oo 



(8) 

h^Yl hn(x, t, X, T) 0^Y1 ^n(x, t, X, T) 



n=0 n=0 



Solenoidality conditions for the flow, magnetic fleld and vorticity result in the follow- 
ing equations: for any n > 0, 

Vx • Mh = Vx • (h„)^ = 0, (9a) 

Vx ■ {V„};, + Vx ■ {^n-l}h = 0, (9b) 

Vx ■ {K}h + Vx ■ {K-i}h = 0, (9c) 

Vx-{u;„}, + Vx-{u;„-i}. = 0. (9d) 

(In differential operators with the indices x and X differentiation in the respective 
spatial fast and slow variables is performed; Vx = {d/dXi, 8/8X2, 0). Differentiation 
in fast variables only is assumed in the definitions (6) of linearisation operators. All 
quantities for n < are zero by definition.) 

The identity (7f) translates into the equation 

Vx X v„ = a;„ - Vx X v„_i. (10) 

Equations (9b) and (9c), together with (10) where the index is n changed to n — 1, 
implies that the right-hand side of (10) is solenoidal in fast variables. The solvability 
condition for (10) is 

{U}n)v = Vx X (V„_i)fe. (11) 

This equation has a solution, if the mean of {u}n)v over the plane of slow variables 
vanishes. Equation (10) then reduces to 

Vx X {Vn}h = {i^n}v - Vx X {Vn-l}h, (12) 

and a solution {v„}/i can be found as a sum of solutions to two Neumann problems 
(see Zh08). 

A hierarchy of equations (A.1)-(A.3) (derived in Zh08 and presented here in Ap- 
pendix A for reader's convenience) emerges at different orders of e upon substitution 
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of the series (8) into (7a)-(7c) and expansion of the resulting equation in power series 
in e. Wc will solve them successively, essentially following the procediu^e described in 
Sections 5-9 of Zh08. For n > the order of calculations is not quite straightforward, 
for each n it involves the following steps: 

1°. Derive a PDE in slow variables in ^u>n-i}v and find {vn)h — i^n}h using the mean 
and oscillatory, respectively, part of the spatial average of the vertical component of 
(A.l) at order 

2°. Derive a PDE in slow variables in ((h„_2))?i and calculate {hn)h — i^n}h using 
the mean and oscillatory, respectively, part of the spatial average of the horizontal 
component of (A. 2) at order e". 

3°. Find the potential parts (in fast variables) of a>„, v„ and h„, applying (9). 

4°. Applying results of the previous steps, derive from (A.1)-(A.3) a system of PDE's 
in fast variables in solenoidal parts of {v„}/i, h„ and On- 

5°. Applying two solvability conditions, calculate ((v„)). (The averages considered at 
steps 1° and 2° also represent solvability conditions for the system (A.1)-(A.3), see 
Section 5.) Apply the remaining solvability conditions, if any are not yet satisfied. 
6°. Use the system of PDE's derived at the previous step to calculate a;„, v„, h„ and ^„ 
in the terms of solutions to auxiliary problems and the yet unknown spatio-temporal 
averages of the terms of expansions (8). 

At this stage solution of the system (A.1)-(A.3) for the considered index n is com- 
pleted, and we continue by solving the equations at the next order e'"'^^ starting anew 
with step 1°. Since the equations for n = studied in the next section are homoge- 
nous, they represent a special case, in which the flow of operations is slightly altered. 
In what follows, we indicate at each point, which step is being carried out. Although 
we stop as soon as derive the amplitude equations, further application of the proce- 
dure would yield solutions to an arbitrary number of successive systems (A.1)-(A.3) 
in the hierarchy. Thus a complete asymptotic expansion (8) can be constructed. 

4. Order equations 

Equations (A.1)-(A.3), (10) and solenoidality conditions (9b)-(9d) reduce for n = 
to 



Differentiation of (la), (lb) and (3) in Xk, k — 1,2, or in time t demonstrates that 
S;-'' = {dn/dxk,dU/dxk,de/dxk) and S*'* = {dn/dt,dH/dt,de/dt) are solutions 
to (13). They satisfy the boundary conditions, whichever of them (no stress or no 
slip boundaries, built of perfectly conducting or a dielectric material, isothermal or 
transferring a prescribed heat flow) are chosen for the flow, magnetic fleld and tem- 
perature. This is a consequence of the assumed spatial and temporal invariance of 
the equation, governing the CHM regime V, H, O. The three solutions to (13) are 
not guaranteed to be linearly independent; for instance, S*'* = for steady CHM 
regimes. Also, S^'"^ are linearly dependent, if the CHM state is independent of a 
horizontal variable aiXi + a2X2 for some constant ai and a2, not vanishing together; 
we will not consider such essentially two-dimensional CHM regimes. 



u>o = Vx X Vo, 



£(u>o, Vo, ho, 6*0) = 0, 

V • Wo = V • Vo = V • ho = 0. 



(13a) 
(13b) 
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Whether other solutions to (13) exist depends on the assumed boundary condi- 
tions and on the parameter values. In what follows wc consider boundary conditions 
(2a), (2c) and (2d). They are inherently consistent with the solenoidality conditions 
and thereby they are widely used in numerical studies of thermal convection and 
convective dynamos (see, e.g., Podvigina 2006, 2008b). 

Following Zh08, we introduce the operator of linearisation in the form 

M{u:,h,e) =C{u},nuj,h,e). 

Here 71 is the inverse curl, acting from the space of solenoidal vector fields globally 
bounded in the layer and satisfying (2b), with a zero spatial mean of the vertical 
component, into the space of solenoidal vector fields satisfying (2a), with a zero 
spatial mean of the horizontal component (see ZhOS, Section 3). The domain of Ai, 
which we denote by T>, consists of vector fields (a;, h, 9), globally bounded in the layer, 
satisfying the respective boundary conditions, and such that uf and h are solenoidal 
with (uj)^ = 0. 

The operator adjoint to M is M* = {{M*)'^, {M*)^, {M*y), 

{MTi^, ^' ^) = ^ + ""^'^ + 7^'{-V X (V X (V X w)) 

+H X (V X h) + ^^ X (V X a;) - r— + ^^63 - ^VG}/,, 

0x2. 

dh. 

{M*f{<jj,\i) = — +r;V2h+Vx(Hx(Vxa;))+P((Vxu;)x(VxH)-Vx(Vxh)), 

(M*)^(uj, h,e) = — + KV^e + (V • V)^ + (3e3 • (V X w). 

Here the operator TZ' is another "inverse curl" : TZ'v = a? is a solution to the problem 

V^w = -V X V, V • uj = 0, (u;)^ = 

satisfying (2b), and V is the projection of a three-dimensional vector field into the 
subspacc of solenoidal fields with a zero spatial mean of the horizontal component, 
satisfying boundary conditions (2c) (see Section 3 in ZhOS for a detailed explanation). 
The operators Ai and Ai* turn out to have the same domain. 

Clearly, A1*(0, C^,0) = 0, where is a constant horizontal vector. Conse- 
quently, there exist eigenmodes S* = (S'^, S'*, S^) of the operator M. associated with 
the zero eigenvalue, for which ((S''))/i 7^ 0. Since spatio-temporal means of S*,'"^ and S*'* 
are zero, the latter neutral modes are distinct. We denote by S^''*^ = (S|^'^, S^!''^, S^'^) 
the neutral modes with non-zero horizontal magnetic components, normalised so that 

Step 2° for n — 0. Averaging in space the horizontal magnetic component of (13a) 
demonstrates that the spatial means of the horizontal magnetic field components 
are time- independent; thus (ho)h = {ho}h- The normalisation is equivalent to the 
conditions (S^''*)/i = at any time t > 0. Spatial means of horizontal components 
of the fiows S^''' = TZS'^'^ are zero. 
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Thus, we have identified up to 5 different solutions to (13) in V existing for 
any parameter values: = S^'^ {k = 1,2), S^_,_2 = S^''^, and if the CHM regime 
V, H, O is unsteady S5 = S*'*. Consequently, the kernel of Ai* is spanned by 4 or 5 
eigenfunctions S^*: = (0, e^, 0) for /c = 1, 2, and S5* exists for non-steady CHM 
regimes. S^* can be chosen to be biorthogonal to S* and normalised so that 

((S^ • S-)) = 5^ (14) 

for 1 < k,j < K. Here K is the dimension of ker A^, Sj is the Kronecker symbol, and 
the scalar product of 7-dimensional vector fields (S'^, S'^, S^) is assumed. The means 
(SDft vanish for all A; < 5, (Sfjh = for k = 1, 2, 5. 

Step 6° for n — 0. By hnearity a solution to (13) can be expressed as 

K 

(u;o,ho,^o) = £S + Eco,feSfe (15) 

k=l 

(interpreted as an equality of 7-dimcnsional vectors). Here ^o(^) -^j ^) — (^O' ^O' ^0) 
is a transient also satisfying the equations (13): initial conditions for ^q, which can be 
found from (15) at t = 0, must belong to the stable manifold of the perturbed CHM 
state V, H, 0, i.e. must exponentially decay in time t. Our task is to construct a 
closed system of equations for the leading order amplitudes co,fc(X, T"). 

If the CHM regime V, H, 6 possesses a symmetry about the vertical axis or parity 
invariance, perhaps involving a time shift T > (see the definitions in Section 7 in 
Zh08), then symmetric and antisymmetric sets of fields are invariant subspaces for the 
linearisation operator A4. Thus any eigenmode is a symmetric or antisymmetric set of 
vector fields (if exponentially decaying transients of arbitrary symmetry are ignored 
in the case of time-dependent modes). In particular, all Sj^ are antisymmetric sets, 
except for S*'*, which is symmetric. Consequently, one can choose a basis in keiAi*, 
where all vector fields are antisymmetric sets, except for a symmetric set S|^, if the 
CHM state V, H, 6 is unsteady. 

5. Solvability conditions 

In this section we consider the system of equations 

M{u:, h, 9) = (P, f^ f), V, • a; = Vx • h = V, • P = V, • = 0, (16a) 

(^). = 0. (16b) 
A solution to (16) exists only, if the following solvability conditions are satisfied: 

= (17) 

(stemming from (16b) ), 

rK = (18) 

and 

((f • SI')) = (19) 
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for A; = 1, 2, 5. Equation (18) is identical to relations (19) for A; = 3, 4. (If the CHM 
state V, H, © is steady, then (19) is not considered for k — 5.) 

We assume henceforth that the invariant subspace spanned by S^, I < k < K, 
is the complete kernel of Ai, where K is the number of eigenmodes in the basis in 
ker A1: X = 4 or 5 for the steady and non-steady CHM states V, H, 6, respectively. 
In the absence of space periodicity the conditions (17)-(19) may be insufficient for 
existence of a globally bounded solution to (16): already quasi-periodicity in hor- 
izontal directions can be problematic. Nevertheless we do not impose periodicity 
conditions, which is too restrictive, but assume that the auxiliary problems of the 
form (16) do have solutions in "D, if (17)-(19) hold true. Subtracting an appropriate 
hnear combination of S^, A; = 3, 4, which belong to ker A4, one obtains a solution to 
(16) with ((hK = 0. 

For n = 0, (11) reduces to {u}o)v — 0, which is compatible with (13) by virtue of 
the identity 

(£-(a;,v,h,^)), = -^. (20) 

The solvability conditions can be used to demonstrate that, unlike in the case of 
forced convection, (vo)/i = 0. 

Step 1° for n = 0. From (A. 4) for n = 1, 

and hence (11) for n = 1 together with (9a) for n = implies 

However, ( {vo}/i|x=£(a;i,a;2)) is asymptotically smaller than any power of e (see Zh08, 
Appendix C). Thus (7g) requires Vq — 0, i.e. 

(vo)h = ((voK, (21) 

and vanishing of the average of ((vo))^ over slow spatial variables. 
Step 3° for n = 0. Consequently, (13a) reduces to 

M{uJo, ho, ^o) = (((vo)),. • V^){n, H, 0). (22) 

Step 4° for n = 0. Scalar multiplying this equation by S** and averaging the re- 
sult, we find that 

((sr-((((vo)).-Vx)(i^,H,e)))) = o 

is a necessary condition for existence of its solution. Hence, biorthogonality (14) of 
S;^ to S** for k,j = 1,2 implies, together with (21) and (11) for n = 1, 

(vo)ft = (wi), = 0. (23) 
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In other words, solvability of (22) in V requires that the mean perturbation flow 
vanishes. We further discuss this issue in Section 10. 

6. Order and £^ equations: a— effect in the leading order 

In this section we derive the a— effect operators, present in the CHM system under 
consideration. 

Step 1° for n — 1. Substitution of the flow and magnetic field (15) into (A.4) for 
n = 2 yields 



Vxx(EE<.^ + rl, (24) 

where 



= {SlH^ + n{s^^u - SlVrn - Y{SIUh, (25) 
I- = (V X (Vx X il) - VVx ■ - H X (Vx X to) + HVx ■ t.W 

The differential operator in the right-hand side of (24) represents the AKA-efFect 
(anisotropic kinematic a— effect) operator in the leading order. Integrating (24) in 
fast time, find 

= +v^xij:j:{{[ dt}^+ {{ [ cdt }] . (26) 



. k=l m=l 







Thus {u}2)v is well-defined, if the means ((/q aj^ ^ dt)) are. This is the condition of in- 
significance of the AKA-effect in the leading order. It is stronger than the solvability 
condition ((f'^))t, = for the equation C'^{u),v,h.,9) — which follows from (20) 
and in the case of (A. 1) for n = 2 is 

«J = 0. (27) 
The solution to equations (11) for n = 2 and (9a) for n = 1 is therefore 

where Vx^ is the inverse Laplacian in slow variables, acting from the space of globally 
bounded functions, whose average in X is zero, into the same space; 

r = {{ ^'dt }} , ^'(x, t, T)^e- VxVi^(Vx • r ). 

Step 2° for n = 1. A possible presence of magnetic a— effect in the CHM system 
with electrically conducting boundaries (2c) can be revealed similarly. Averaging of 
the horizontal component of (A. 2) for n = 1 in fast spatial variables with the use of 
the identity 

(£V,h)). = -^ (29) 
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and substitution of the flow and magnetic component of (15) yields 



^-Vxx(f:a^co,. + |'^), (30) 

where 

= (V X + S", X H)„ (31) 

I'^^d^xH + Vx^,. 

The differential operator in the right-hand side of (30) represents the magnetic 
a— effect. The equation yields 



(hi). = ((hi)), + Vx X f f: {{ /J cxl dt }} CO,, + {{ [~e<it }} ) . (32) 



Magnetic cy— effect is insignificant in the leading order, if the means ((/q a^dt)) ex- 
ist and as a result (hi)/j is well-defined by (32). The condition of insignificance of 
the magnetic a— effect is stronger than the solvabihty condition ((f'*))/t = for the 
equation C^{y, h, 9) — f'*, which follows from (29) and in the case of (A. 2) for n = 1 
is 

io^D) = 0. (33) 

If the perturbed CHM state is steady or periodic in time, the 4A' scalar relations 
(27) arc sufficient for insignificance of kinematic a— effect, and the K scalar relations 
(33) are sufficient for insignificance of the magnetic a— effect. Elementary algebra 
shows, that these conditions hold true for A; = 1, 2, 5: 

= = for ^ = 1, 2, = ^ (H//„ - VK^),, = ^ (V x H),. 

If the CHM regime V, H, O has a symmetry considered in Zh08, then (27) and (33) 
are satisfied for all k, and thus the AKA- and magnetic a— effects are insignificant in 
the leading order. If the symmetry is spatial (i.e. the time shift T is zero), then a* = 
for aW. k < K = A\ if it is spatio-temporal, then for all k < A OL\(t + T) = — Q:*(t) 
and ^/ga'dt]} satisfies the same relation (see the proof in Zh08, Section 7); the 
remaining vector fields ol\ for A; = 5 satisfy OL\{t + T) = cx\{t), and hence -^/q cx\ dt^ 
also satisfy the latter relation. 

The terms f/o I'^dil in (26) and f Jo |'*di> in (32)_are not problematic, because 
^'^ and exponentially decay and thus the means ((/q ^'^dt)) and ((/q ^'*dt)) are well- 
defined. The quantities ■8^/o?'^di^, and ^fok^'^t^ also decay exponentially (see 
Zh08, Appendix D). 

The multiscale approach remains feasible even, if (27) and/or (33) do not hold 
true. It is well-known (see DubruUe and Prisch 1991), that in this case another slow 
time scale is appropriate, T ~ at. Then the new term d{ho)h/dT emerging in the left- 
hand side of (30) balances the magnetic a— effect term. Then equations (24) and (30), 
supplemented by (38) (see the next section) if the CHM regime V, H, is unsteady, 
constitute a closed system of linear first-order PDE's (together with the solenoidality 
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conditions (9a), forn = for the mean magnetic field, and for n = 1 for the fiow (37) 
also derived in the next section; in view of (23). (38) and (37) arc not affected by 
the change in the time scaling). This system turns out to be mixed: while (30) is an 
evolutionary equation, (24) and (38) are not. Solutions to such systems generically 
exhibit unbounded exponential growth. (Multiscale expansion can yield amplitude 
equations, where the a— effect coexists together with the molecular diffusivity and the 
nonlinear advection, but this requires construction of the asymptotic expansion along 
a different line in the parameter space. For instance, Frisch et al., 1987, obtained such 
mean-flow equations for perturbations of flows in a hydrodynamic weakly nonlinear 
stability problem in the limit of small Reynolds numbers R, assuming that the scale 
ratios were linked to R.) Consequently, from now on we focus on a potentially more 
interesting case of insignificant ct— effect, where possible growth of perturbations may 
saturate due to nonlinearity. 

7. Order equations: a— effect, insignificant in the leading order 

Assuming that the AKA- and magnetic a— effects are insignificant, in this section we 
solve equations (A.1)-(A.3) for n = 1. 

Step 3° for n = 1. To apply the solvability conditions (17)-(19), this system must 
be transformed to take the form of the problem (16). Thus, the gradient parts of 
a;i,Vi and hi, and the mean part of Vi must be isolated. These quantities can be 
found using (12), (9b)-(9d) for n = 1, (23) and (15): 

K 2 

^i = EJ:7^^mXSl + u:[, (34a) 



k=l m- 



K 2 f) 
k=l m=l ^^rn 

hi = EE^v.4, + h;, (34c) 



k=l ra- 



where 

Vx • = Vx • vi = Vx ■ K = 0, (o^'i). = {v[)h = 0, Vx X v'^ = (35) 
and s^j.(x, i) are globally bounded solutions to Neumann problems 



^x^m,k~ y^kJm, 

'^^m,k 



"5m,fe — {^k}m^ 



= 0; 

X3=±L/2 

= 0. 

X3=±L/2 



Consequently, u>'i,v'i and h^ satisfy the same boundary conditions, as Wi, Vi and hi, 
respectively. 
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Step 4° for n = 1. Upon substitution of (34), equations (A.1)-(A.3) forn = 1 take 
the form 

fc=l m=l C/^m 

- 2i/(Vx ■ Vx)u;o - Vx X (V X u?o + vo X ^^ - H X (Vx X - ho X (Vx x H)) 
-Vx X (-H X (Vx X ho) + Vq X uJo - ho X (Vx x ho)) - /9Vx6'o x 63, (36a) 



CCo^A 



M\u;'„ h[) = - E E ^'(Vx<„ Vx4.) ^ + ((vi), ■ Vx)H 

fe=l m=l 



-277(Vx • Vx)ho - Vx X (vo x H + V X ho) - Vx x (vo x ho), (36b) 

K 2 Q 

M\u;[, ^1) = - E E >C'(Vx<„ 0) ^ + ((vi). • Vx)e 

k=l m=l iJ^m 

-2«(Vx • Vx)^o + (V • Vx)^o + (vo • Vx)^o. (36c) 

Solvability condition (17) is clearly verified for (36), and can be determined 
from (35). Condition (18) yields expression (32) for (h.i)h (sec Section 6). 

Step 5° for n = 1. Scalar multiplying (36) by S**, j = 1,2, using their biorthogo- 
nality to Si and normalisation (14), and employing (15), (28) and the identity 

-^+^V^S^+VxS-+S",xfi-Hx(VxS^)-S^x(VxH)+rS^xe3+/35^3 = V5,^ 
obtained by "uncurling" the vorticity component of (13a), we find 

(37) 

where it is denoted 

/3.^(((E:-sr)),((E:.sr)),o), 

and E* are 7-dimensional vector fields: 

E-,,^,, = £(0, VxC,., Vx«:^,„ 0) - ({{ 1^ < dt } ■ Vx) (n, H, e) 

+ Vx X (V X (e^ X S^) - H X (e^ x S^)) + x Vx^^ + r{Sl)^es, 

f)Gh ace \ 

277^ + e^x(VxS^ + S^xH), 2«:-^ - K.^^ , 

^'2,m,k ^ Pm,k (Vx X (S", X S^^ - X (Vx X S^) X - X (Vx X S^)), 

Vx X {s>l. X + X S^j , — (S^ ■ Vx)>S'^ — (S^ ■ Vx)*?^) , 
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^h,,rn,k = P^, (({{ icy^m,k) , }} e, + {{ «,)^ dt }} 6,) ■ Vx) (^^, H, 6), 

p^ fc = 1 for m < /c, p„, = 1/2 for m = k, and Pm,A; = for m > k. If the perturbed 
CHM regime V, H, G possesses a symmetry considered in Zh08, vector fields E* for 
k < A are symmetric, Sj^' for k — 1,2 are antisymmetric, and thereby {vi}h = 0, if 
V, H, © is a steady CHM state. 

If the CHM regime V, H, G is unsteady, scalar multiplication of (36) by Sg*, with 
the use of its orthogonality to S^, k — 1, 2 and normalisation (14), yields 

(38) 

where 

/3: = ((E:-sr)). 

All the solvability conditions for (36) have been satisfied, and thus a solution can be 

constructed. 

Step 6° for n = 1. Let V{sl'^ ,aL^,a^) denote the projection of a vector field 
(a"^, a'*, a^) to the subspace, orthogonal to ker Jii* (in the sense of the scalar product 
(((a^, a^, af) ■ (a^, a^, a^)} ). If the CHM state V, H, © is unsteady and symmetric, 
application of P amounts to projecting out the S5* component, since the right-hand 
sides in (36) are symmetric, and S5* is the only symmetric eigenmode in kcrA^*. 
In view of relations (15) and (28), and due to linearity of this problem, it has the 
following solution: 



k=l \ 



+ J2 f ^m,k-^^ + ^lj,m,kr.y JTyTy I ^ Q'm,kC0,kC0,m ) ■ (39) 

Here G;^^^ = (G;:;,^^, G"^^^, Gl^^^, G^^fe) solve auxiliary problems of type II : 

M{G'^^^) = V^l^^,, (40a) 

V. X Gl,^, = G;;,,, (40b) 

V. • G- = V. • g;;,,, = V. • Gt^, = o. (40c) 

Vector fields Q;, ^ = (Qm,fe> Qm,fc> Qm,fe> Qtn,k) solve auxiliary problems of type III : 

M-{Q'^,,) = VEl^^,„ (41a) 

Vx X Q;;,, = Q:;,,, (41b) 

V. • q:;,, = V. • q:;,, = Vx • q::,,, = o (4ic) 

(Qm,fc = for m > /c). 
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Vector fields ^ ^ = {^tj,m,k^ ^Um^k^ ^tj,m,k^ ^f,j,m,k) solve auxiliary problems 
of type IV: 

M ( YV,^,,) = PE-,,,^. (42a) 
Vx X Y^,,,^,, = Y^,^^^„ (42b) 

(Y*j,™,fc = Ofori> j)- 

Finally, |* = (^^ , H) solve the problem 

>C"(ir, li, e?) = -2^(Vx • Vx)eo - Vx X (V X + 1^ X 

-H X (Vx X to) -t^x (Vx X H)) - Vx X (vo x + x (a;o - C) 

-H X (Vx X - ho X (Vx X t^) -ilx (Vx X (ho - il))) - /^VxCo x eg, (43a) 
>C'(£L C?) = -2r?(Vx • Vx)^? - Vx X (^^ X H + V X ^o') 

-VxX(voX^^' + |[;x(ho-|o')), (43b) 

il) = -2«:(Vx • Vx)eo + (V ■ Vx)eo + (vo ■ Vx)eo + (Co " Vx)(^o - Co'), (43c) 

VxX^^-|r = -Vxx|^, (43d) 

Vx-Cr = -Vx-£^, Vx-|^ = -Vx-C Vx • = -Vx • Co'- (43e) 



^'m,k^Qm,kJ^i,j,m,k C* must Satisfy the boundary conditions (2a)-(2c) and 
(4). Spatial means of their horizontal flow components and of vertical vorticity 
components must vanish, as well as spatio-temporal means of horizontal magnetic 
components. 

Together, (40c), (41c), (42c) and (43e) are equivalent to (9b)-(9d) for n = 1. It 
suffices to impose these conditions for the vorticity and magnetic parts of a solution 
at t = (this can be deduced taking the divergence of the vorticity and magnetic 
equations of the problems (40)-(43) ). Equations (40b), (41b), (42b) and (43d) are 
equivalent to (12) for n = 1. 

The spatial means of the initial conditions for these problems can be determined 
integrating the equations for vorticity and magnetic perturbations and considering 
the spatio-temporal means of the results: 

{G'k)H^_=-(emX tcx^dt)), (44) 



t=o W ™ Jo " II 

t _ 



(the means exist by the assumption that the magnetic a— effect is insignificant); at 
any t > 

{^m,k)v — {Qm,k)v — 0^i,j,m,k)v — {^l)v — {Qm,k)h — 0^i,j,m,k)h — 0- (45) 
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Averaging the magnetic part of (39), we obtain 



t=0 



t=0 

k=l m=l 



T=0 



Now initial conditions for the problem (43) can be determined from (39) and (34) at 
t = 0. 

The choice of initial data must ensure that all solutions to the auxiliary problems 
are globally bounded with their derivatives. Any solution is globally bounded for any 
smooth initial conditions, if the perturbed CHM state V, H, is space-periodic and 
linearly stable to small-scale perturbations (sec Appendix B in Zh08). Modification 
of initial conditions for G^^, Q^ ^ and Y*^^^ into other permissible ones implies 
the respective changes in initial conditions for These changes must belong to the 
stable manifold of the perturbed CHM state V, H, ©, so that the resultant changes 
in solutions to the auxiliary problems and in ^* decay exponentially in fast time. 

Since decay exponentially in fast time together with derivatives (see Section 4), 
the right-hand sides of equations (43) also do. If the CHM state V, H, 9 is hnearly 
stable to small-scale perturbations, this implies that ^' decay exponentially, and any 
changes in j,, Q^,fc and Y* j „ due to a permissible variation of the initial data for 

also exponentially decay (see Zh08, Appendix B). Otherwise, exponential decay 
of ^* must be ensured by an appropriate choice of the initial data. 

If the CHM regime V, H, 6 is steady or periodic in time and periodic in spatial 
variables, generically one can find Sj., G^ j,, ^ and Y*^-^^, which are steady or 
have the same periodicity, respectively (sec Section 4 in Zh08). If the perturbed 
CHM regime has a symmetry of the kind considered in ZhOS, arc antisymmetric 
sets, except for Sj, which is a symmetric set of vector fields for non-steady regimes. 
Consequently, G^ ^ and Y'^^^^ are then symmetric sets for A; < 5 (by construc- 
tion the antisymmetric part of any permissible solution to the problems (40)-(42) for 
A; < 5 exponentially decays and hence is irrelevant), as well as Q'^^^; G^ g and Y'^^ g 
are antisymmetric sets. For symmetric sets conditions (44) and (45) are automati- 
cally satisfied (except for the mean of the vertical component of vorticity, vanishing 
of which is not implied by parity invariance of the perturbed state) . If the perturbed 
CHM state is steady or its symmetry is without a time shift, the right-hand sides of 
auxiliary problems (42) are zero and hence Y^'^^ = 0. 

8. Order £^ and equations: amplitude equations 

Step 1° for n — 2. Combining the flow velocity component of (39), (37), (34b) and 
(28), we flnd 

+ E Qm,kCo,kCo,rr}\ , (46a) 

, m=l / 

where ^ 

^m,k — ^m,k + (^l,m,k + "fi / "m,fe J + ^^^m,k^ 
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,m,ki 

Similarly, 

^ ^ £ (^''''' ^ £i {^'^''^ ^ % S al^^, 

+ E Qm,feC0,fcC0,^) , (46b) 
m=l / 

where 

T^/i _ r^h _|_ VT 

Averaging (A. 4) for n = 3 in fast time, substituting (23), (27), (46) and the flow 
and magnetic components of (15), and recalling that Ci decay exponentially, 

we obtain 

^ ^ d ( ^ d ( V ^ ^ V 9'^ -2 ^ 

X ^ ^ I \T)l „^ i.Co^k + E E ^l,i,j,m,k QX dX^^^^'^ 

k \ 

+ E K,m,kC0,mC0,k] =0, (47) 
m=l / 

where ^ ^ 

^n,m,k = (("^GJ^^j^. — V(G^^^)^ + HnG^f, + H(G^^^)^))/i, 
*^n,ij,m,fc ^ ((~^^i,j,m,fe ~ ^(■^j,m,fe)„ "I" ^ri^ i,j,m,k "I" -^(■^j,m,fe)„))'i' 

Step 2° for n = 2. Averaging the horizontal component of (A. 2) for n = 2 in fast 
variables, substituting (46) and flow and magnetic components of (15), recalling that 

and ^* decay exponentially, and taking into account relations (29), (33) and the 
boundary conditions for V, H, v„ and h„, we find 

-^((hoK + r;V^((hoK + Vx X E ( E ^{^l,kCo,k 

k=l \ m=l ^ \ 
2 j Q2 ^ \ k \ 



(48) 



where 



Ka,k = ((V X q:^,, - h X q;;,, + s^, X s:^)),. 
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D* are coefficients of tfie second-order operators representing tfie so-called aniso- 
tropic combined eddy dijfusivity correction, d* are coefficients of pseudodifferential 
operators, formally also of the second order, which can be regarded as representing an 
unconventional non-local anisotropic combined eddy diffusivity correction. All d* = 0, 
if the perturbed CHM state is steady or possesses a symmetry considered in Zh08 
without a time shift. A* are coefficients of quadratic terms representing the so-called 
combined eddy advection correction. 

Equations (47) and (48) are solvability conditions for the systems (A.1)-(A.3) 
for n = 2 and 3. Together with (38), if the CHM regime V, H, 6 is unsteady, they 
constitute a closed system of equations for the leading terms of mean perturbations 
(note that Co,jfc+2 = ((ho))jfc, k — 1,2, since (S^^2)/i — ^fe)- The vertical component of 
(48) and the horizontal component of (47) vanish identically (i.e., (48) represents two 
scalar equations and (47) just one). 

Unlike for forced CHM systems, the system of amplitude equations is mixed, in- 
volving both evolutionary and non-evolutionary equations in slow variables. Equation 
(48) for the mean magnetic field perturbation is evolutionary. It preserves solenoidal- 
ity of the mean magnetic perturbation (condition (9a) for n = 0) in slow spatial 
variables, which therefore becomes just a constraint for the initial condition. Equa- 
tion (47) resulting from the vorticity equation does not bear any similarity with the 
original equation - both the derivative in slow time and molecular diffusivity operator 
are absent. Equation (38), emerging if the CHM regime V, H, © is unsteady, is also 
non-evolutionary, as well as the condition of solenoidality in slow variables (9a) for 
the mean flow ((vi))/j (37). The system of amplitude equations is still underdeter- 
mined, if the CHM state is steady and possesses the symmetries considered in Zh08, 
resulting in ((vi))/i = 0. The missing equation is then the condition of solenoidality 
in slow variables of the mean fiow {v2)h- 

9. Mean flow perturbation {v2)h for a symmetric steady CHM state 

In this section we calculate {'V2)h ioi a steady CHM state V, H, G, which is parity- 
invariant or symmetric about a vertical axis. In this case, in view of (23), (26), (28), 
(37) and (11) for n = and 1, 

(vo)fe = (vi)fe = (a;o)^ = {uji)v = {^2)v = 0. 

This is not an independent calculation - we advance by two steps in sohition of the 
system (A.1)-(A.3) for n = 2, following the general procedure for treatment of the 
systems in the hierarchy that is outlined in the end of Section 3 (although the solution 
itself is not required for our purposes). 

Step 3° for n = 2. We need to apply solvabihty conditions (19) to equations (A.l)- 
(A.3) for n = 2. To do this, wc transform the system to fit the form of the problem 
(16) by isolating the gradient parts of W2, V2 and h2, and the mean part of V2 using 
(12), (9b)-(9d) forn = 2 and (39): 

X Q;;,;. + u;^, (49a) 

fe=l n=l '^^n \ m=l rn=l / 
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K 2 / 2 K \ 

fe=ln=l "^^^ \m=l '^^"i m=l / 

K 2 / 2 K \ 

h2 = E E ^ E E co.^co,.V.<„,, +h^, (49c) 

k=l n=l \ m=l rn=l J 

where 

Vx • ^2 = Vx • = Vx • h'2 - 0, (u;^), = (V^);, = 0, Vx X = J^, 
gi* „j j(j(x, t) and g* „j j.(x, t) are globally bounded solutions to Neumann problems 



^x9n,m,k ~ ~i^m,k)n: 
~{Qm,k)ni 



dxs 



m.k 



rln,m,k 



,m,k 



8X2, 



X3=±L/2 



X3=±L/2 



X3=±L/2 



^x1n,m,k^ {Qm,k)ni 



m,k 



8X3 



0; 

0; 
0; 
0, 



X3=±L/2 

and in agreement with the previously used notation (recall that /3, = in (37) for a 
perturbation of a symmetric CHM steady state) 

^m,k — {^m,k + X S^, G^^^ + Vx^^^fc, G'^ ,^ + VxS^^fc, Gm,k)- 

The fields a;2, V2 and satisfy the boundary conditions for vorticity, flow and mag- 
netic field, respectively. 

Step 4° forn = 2. Substituting (49), (15) where Co = 0, and 



K 



da 



(u;i, vi, hi, ^1) = E ( S^ci,, + E Gl^,^ + Y: Q'm,kCo,kCo,. 



K 



k=l 



m=l 



8X, 



m rn=l 



(a consequence of (39) and (34), where = 0) into equations (A.1)-(A.3) for 

n — 2, noting that Ci^k enter these equations only via terms which are symmetric 
sets, scalar multiplying the resultant equations by S**, j = 1,2 and using their 
biorthogonality to Si and normalisation (14), we find 



K / 2 2 



{V2)h = E I E E /^l,n,m,fco^ 2,n,m,k^0,r 



k=l \m=ln=l ^^m^^n ^=1^=1 

K K \ 

+ EE/3 

,n I • 

m=l n=l / 



8X^ 



(50) 
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Here it is denoted 

f^i,n,m,k — i{^i,n,m,k ' ^1 ) i {^i,n,m,k ' ^2 )) 0), 

where E* are 7-dimensional vector fields: 

+ e„ X (V X G-,, + Gl^, X 11 - H X (V. X Gt^,) - G^^, x (V. x H) 
-H X (e„ X S^) - PGl^,es) - V. x (H x (e„ x G^^^J), 

^"^^^'-'^ + ^" ^^^'^ X H + V X G':^^,), ^f^^flk - VnGl^j ; 

E* = p* -I- p* 

2,n,m,A; n,m,k ' n,k,m, 

+ (Vx X (G^j X + X G'^^j^ — G^;, X (Vx x S'^) — x (Vx x G^^^)), 
Vx X (G^^;;. X + X G^^;;.), -{Gl^^ ■ V^)S^ - (S^ • Vx)G^_fe) ; 

V V h f ^ u> 

^n,m,k = ^{^n X Q^_fe, Vx Qn,m,k^ Qn,m,k^ 0) + I '^^'q^^^^^ 

+ e„ X (V X Q-,, + Q;;,, X O - H X (Vx X Q::,^,) - Q^;,,, X (Vx X H) 
+ S;; X S;^ - S::, X (Vx X S^) - (5Ql^,e,) - Vx x (H x (e„ x Q^^,,,)), 

^^^^^t,k + en X (Q;;,,, X H + V X (tm,k + X S,^), 2«: Ag^^^ _ y„g^^^ j ; 



%n,m,k ^ (Vx X (q;;,,, X s:^ + s:; X - q^;,,, x (Vx x s^^) - s^: x (Vx x q;^,,)), 

Vx X (Q^,A: X S(j + X Qm^fc), — (Qm^fc " Vx)5'^ — {^^n ' Vx)Qm,fc) • 

Thus, for a steady CHM state possessing a symmetry which guarantees insignificance 
of the a— effect, the missing amplitude equation - the solenoidality condition for the 
flow (50) - has been expressed in the terms of the amplitudes Co,fe. 

10. Non-zero initial mean flow perturbation ((vo))/i 

It has been found in Section 5 that unless the horizontal part of the mean flow 
perturbation vanishes, no solution of order £° equations (22) exists in the functional 
space T> (introduced in Section 4). We discuss here the possibility of relaxing the 
condition that the solution belongs to V. 

Direct substitution demonstrates that (22) has a solution 
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It satisfies all the conditions for vector fields in D, except for it is not uniformly 
bounded in the fast time t. (Other solutions to (22) can be constructed, allowing 
their linear growth in horizontal directions, but they seem too unphysical.) 

If growing in time solutions are regarded as acceptable, further application of the 
asymptotic techniques does not require considering dependence of the perturbation 
(8) on the slow time T. Construction of an expansion of large-scale solutions to the 
problem (1) in power series becomes more straightforward. At each order one 
proceeds in two steps: 

1. Calculate {ujn+i)v and {vn)h using (A. 4) at order e""*"-*-. 

2. Find {u)n}v, {'^n}h, hn and 9n integrating in time t the system (A.1)-(A.3) at order 

For the sake of simplicity, suppose that the CHM state fl, H, T is steady and 
periodic in horizontal directions. Then the system (A.1)-(A.3) can be solved ex- 
panding the unknown fields a;„, h„, 9n in a series in eigenfunctions of the operator of 
linearisation of the original equations, £, which is an elliptic operator (not involving 
the derivative in time). If the steady CHM state is stable to short-scale perturba- 
tions, all eigenvalues of this operator have non-positive real parts. The coefficients 
of the expansion are accessed by scalar multiplication of the system by the biorthog- 
onal eigenfunctions of the adjoint operator. By virtue of (51), the zero-order terms 
ujo, vo, ho, 6*0 grow linearly in time. As a result, u)n, v„, h„, 9n experience a polynomial 
growth in time, the degree of the polynomial being 2n -|- 1. 

Thus, free convection is apparently unstable to large-scale perturbations with a 
non-zero mean flow. However, this is not a genuine instability, but rather variation 
(at a given point in space) of the profile of the perturbed CHM fields f2, H, due 
to advection by the mean fiow {{vojh- Were this mean flow constant, one could rely 
on the Gahlean invariance to achieve ((vo))/i = 0, i.e., to overcome the difficulty by 
considering the problem in the comoving coordinate system. (Note that the solution 
(51) to the order problem (22) represents the second term in the Taylor expansion 
in time of the fields O, H, transported by the fiow ((vo))^.) It is conceivable that a 
similar transformation modifies the present weakly stability problem into one solvable 
in V, but dependence of ((vo));j on slow variables makes the new problem significantly 
harder, than the problem considered here. 

11. Conclusion 

We have investigated, in the Boussinesq approximation, weakly nonlinear stability 
to large-scale perturbations of a convective hydromagnetic regime V, H, T in a hori- 
zontal layer rotating about a vertical axis. We have assumed that the CHM system 
is free, i.e., there are no source terms in the governing equations (1) and bound- 
ary conditions (2). Such system is translation invariant in space and time. We have 
derived expressions for the a— effect tensors (25) and (31). If the CHM regime is sym- 
metric about a vertical axis with a time shift, or parity invariant with a time shift, 
the a— effect tensors do not contribute to the leading order amplitude equations, 
i.e. the a— effect is insignificant in the leading order. Applying the homogenisation 
techniques within the approach developed in Zh08, under the assumption that the 
a— effect is insignificant in the leading order (because of the symmetries or otherwise) 
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we have derived a system of amplitude equations, (47) and (48), augmented by (38) 
if the perturbed CHM regime is non-steady. This system, governing the evolution of 
large-scale perturbations in slow time, is closed, when considered together with the 
condition of solenoidality in slow spatial variables of the leading term of expansion 
of the mean flow perturbation. In a generic a— effect-free setup, this is the condi- 
tion of solenoidahty of the flow (37). If the CHM state is steady and possesses a 
symmetry guaranteeing insignificance of the a— effect, ((vi));^ = and this condition 
becomes trivially satisfied; then the remaining amplitude equation is the solenoidality 
condition (9a) for the flow (50). 

Similarly to the case of forced thermal hydromagnetic convection, which was stud- 
ied in Zh08, the system of amplitude equations involves a linear operator of combined 
eddy diffusivity correction and a quadratic operator of anisotropic combined eddy ad- 
vection. Generalised diffusivity correction is described by non-local pseudodifferential 
operators, formally of the second order, like the standard diffusivity. All the opera- 
tors are anisotropic. No similarity of amplitude equations for free CHM regimes and 
mean-field equations for forced CHM regimes is observed beyond this point: Because 
of the differences in the structure of the kernels of the operators of linearisation in the 
vicinity of the perturbed regime in the two cases, the system of amplitude equations 
that we have obtained here does not involve mean flow perturbation, and it is mixed. 
Whereas equation (48) for the mean magnetic perturbation is evolutionary, hke it is 
for forced convection, the remaining ones, (47) and (38), involve neither time deriva- 
tives, nor molecular diffusivity operators. If the CHM state is steady and possesses a 
symmetry guaranteeing insignificance of the a— effect, the solenoidality condition for 
the fiow (50) reduces to a non-evolutionary third-order partial differential equation 
with a cubic nonlinearity. 

As already mentioned above, the neutral linear stability modes S^'^ and S*'* exist 
for free CHM regimes for any combination of physically sound boundary conditions 
because of the spatial and temporal invariance of the CHM regime V,H, 0. This 
opens a possibility of a similar analysis of linear or weakly nonlinear stability of 
CHM regimes for boundary conditions different from the ones considered here. For 
instance, it is feasible for electrically conducting fiuid contained in a horizontal rotat- 
ing layer confined between half-spaces of dielectric material with isothermal no slip 
boundaries, when the collection of neutral modes is the smallest one - for generic sets 
of parameter values it is comprised of the three above-mentioned modes. Note that a 
direct application of this approach to the study of stability of forced convection fails 
for these boundary conditions. 

It is possible to perform the same analysis for branches of CHM regimes emerging 
from symmetric ones near Hopf and pitchfork bifurcations, like it was done in Zh08 
for forced CHM regimes. New a— effect terms emerge in the amplitude equations 
for the evolution of large-scale perturbations in stability analysis of regimes near the 
point of a symmetry-breaking Hopf bifurcation, and an additional amplitude equation 
involving cubic nonlinearity emerges in stability analysis for regimes emerging in 
a symmetry-breaking pitchfork bifurcation. Since the derivations, although quite 
straightforward, are very technical, we have not carried out them here. 

We have not prescribed any boundary conditions for the amplitude equations. 
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since this would not affect our derivations. The basic requirement is that the pertur- 
bation is globally bounded - this is the logical foTindation for application of asymptotic 
methods employed here. In computations it is natural to assume periodicity of the 
amplitudes in slow variables. Results of a numerical study of large-scale perturbations 
under this condition, for a class of CHM regimes symmetric about the vertical axis 
and stable to short-scale perturbations, will be reported in a sequel to the present 
paper. 

Acknowledgments 

Part of this research has been carried out during my visit to the School of Engineering, 
Computer Science and Mathematics, University of Exeter, UK, in January - April 
2008. I am grateful to the Royal Society and the University of Exeter for their 
financial support. My research visits to Observatoire de la Cote dAzur were financed 
by the Ministry of Education of Prance. I was partially supported by the grants 
BLAN07-2 183172 (research project OTARIE) from Agence nationale de la recherche 
(Prance) and 07-01-922 17-CNRSL^ from the Russian foundation for basic research. 

References 

DubruUe, B. and Prisch, U. Eddy viscosity of parity-invariant flow, Phys. Rev. A, 
1991, 43, 5355-5364. 

Prisch, U., She, Zh.S. and Sulem, P.L. Large-scale flow driven by the anisotropic 
kinetic alpha effect, Physica D, 28, 1987, 382-392. 

Podvigina, O.M. Magnetic held generation by convective flows in a plane layer, Eur. 
Phys. J. B, 50, 2006, 639-652. 

Podvigina, O.M. Instability of flows near the onset of convection in a rotating layer 
with stress-free horizontal boundaries, Geophys. Astrophys. Fluid Dyn., 102, 2008a, 
299-326. 

Podvigina, O.M. Magnetic field generation by convective fiows in a plane layer: the 
dependence on the Prandtl number, Geophys. Astrophys. Fluid Dyn., 102, 2008b, 
409-433. 

Podvigina, O.M. On stability of flows near the onset of convection in a layer with 
stress-free boundaries, Geophys. Astrophys. Fluid Dyn., 2009, submitted. 

Zheligovsky, V.A. On the linear stability of spatially periodic steady magnetohydro- 
dynamic systems with respect to long-period perturbations, Izvestiya, Physics of the 
Solid Earth, 39 (5), 2003, 409-418. 

Zheligovsky, V. Mean-fleld equations for weakly nonlinear two-scale perturbations of 
forced hydromagnetic convection in a rotating layer, Geophys. Astrophys. Fluid Dyn., 
102, 2008, 489-540. 

Zheligovsky, V.A., Podvigina, O.M. and Prisch U. Dynamo effect in parity- invariant 
ffow with large and moderate separation of scales, Geophys. Astrophys. Fluid Dyn., 
95, 2001, 227-268. 



23 



Appendix. The hierarchy of equations for weeikly nonUnear perturbations 

The following equations arise at order after substitution of the series (8) for the 
perturbation into (7a)-(7c) and expansion in power series, for the perturbed CHM 
state ft, V, H, O independent of e: 



9n) - + ly (2(Vx • Vx){a;n-i}. + V^a;,_2 

+Vx X (V X LJn-l + Vn-1 X - H X (Vx X h„_2 + Vx X h„_i) - hn-1 X (Vx X H) 

n-2 

+ ^(Vfe X UJn-2-k - hfc X (Vx X h„_2_fc + Vx X h„_3_fc))) 
fe=0 

n-1 

+Vx X (^(Vfe X UJn-l-k - hfc X (Vx X K-i-k + Vx X h„_2_fe)) 
fe=0 

-H X (Vx X h„_i)) + /3Vx^n-l X eg = 0, (A.l) 

C\v^, K) - + V (2(Vx ■ Vx){h„_i}, + V^h„_2) 

n— 2 n—1 

+Vx X (v„_i X H + V X h„_i + ^ Vfc X h„_2-fc) + Vx x ^ Vfc x K-i-k = 0, (A.2) 

k=0 k=0 

an 

/:'(V„, - + K (2(Vx ■ Vx)^^n-1 + V|^n_2 



n—1 n— 2 

-(V • Vx)^n-l - E(Vfc ■ Vx)^n-l-fe - EK • Vx)^n-2-fc = 0. (A.3) 
fe=0 A;=0 



Averaging the vertical component of (A.l) in fast spatial variables, taking into 
account relations (9b), (9c) and (10) (valid for any index n), (20) and the boundary 
conditions for V, H, Vj and hj, we obtain 

d{uJn)v , d{uJn-2)v ^2 / \ 

^^ + ^^ = -Vx(a;._2). 
+ Vx X (V X (Vx x v„_2) - VVx • v„_2 - H x (Vx X h„_2) + HVx-h„_2 

n-3 

+ J2(^kX ( Vx X v„_3_fe) - Vfe Vx • v„_3_fe - hfe X ( Vx X h„_3_fe) + hfe Vx • h„_3_ 



k=Q 



(A.4) 
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